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ABSTRACT 

We detect a positive angular correlation between bright, high-redshift QSOs and fore- 
ground galaxies. The QSOs are taken from the optically selected LBQS Catalogue, 
while the galaxies are from the APM Survey. The correlation amplitude is about a few 
percent on angular scales of over a degree. It is a function of QSO redshift and apparent 
magnitude, in a way expected from weak lensing, and inconsistent with QSO-galaxy 
correlations being caused by physical associations, or uneven obscuration by Galactic 
dust. The correlations are ascribed to the weak lensing effect of the foreground dark 
matter, which is traced by the APM galaxies. The amplitude of the effect found here is 
compared to the analytical predictions from the literature, and to the predictions of a 
phenomenological model, which is based on the observed counts-in-cells distribution of 
APM galaxies. While the latter agree reasonably well with the analytical predictions 
(namely those of Dolag & Bartelmann 1997, and Sanz et al. 1997), both under-predict 
the observed correlation amplitude on degree angular scales. We consider the possible 
ways to reconcile these observations with theory, and discuss the implications these 
observations have on some aspects of extragalactic astronomy. 
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1 INTRODUCTION 

It is now widely recognized that the sky-projected density of 
high-redshift objects can be altered due to weak lensing by 
intervening mass distribution, on a range of angular scales 
(Narayan 1989, Broadhurst, Taylor & Peacock 1995) The 
effect arises because gravitational lensing distorts the area 
on the sky in the direction of the lens. Behind a lens of a 
positive mass density, the area is stretched out. As a re- 
sult: (1) individual sources subtend a larger area on the sky 
and therefore look brighter due to flux conservation; (2) but 
their number density decreases. If the slope of magnitude 
number counts is steep then the first effect wins over the 
second and the net number density of background sources 
in a flux limited survey is increased; if the slope is shallow, 
the number density is decreased. Thus the amplitude of the 
effect depends on the amount and dumpiness of the lensing 
matter, the relative redshifts of the lenses and sources, and 
the slope of the source number counts in the appropriate 
redshift range. 

The effect of dilution or enhancement of the pro- 
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jected number density is most pronounced for sources with 
magnitude-number counts whose slope deviates strongly 
from a=dlogN /dm of 0.4. In the cases where the slope is 
shallower then 0.4, anti-correlations with foreground lenses 
have been observed or suspected: Rodrigues- Williams & 
Hogan (1994) suggest that the anti-correlation of faint UVX 
objects with clusters (Boyle et al. 1988) are due to weak lens- 
ing rather than dust. Broadhurst (1994) observe a deficit of 
red faint galaxies (a « 0.3) behind the foreground cluster 
Abell 1689. 

When the slope is steeper than 0.4, positive correlations 
are detected. These observations can be roughly categorized 
by the angular scale of correlations. On small scales, 3 — 30", 
there are a number of observations (see Narayan 1991, and 
Hewett, Harding & Webster 1991 for reviews). Though di- 
rect comparison between these is difficult, because each sam- 
ple has its own selection criteria and method of analysis, the 
correlations are thought to be reasonably well understood in 
terms of lensing by individual galaxy dark matter halos. 

On 1 — 15' scales the correlations are due to dark mat- 
ter associated with galaxies, or clusters of galaxies, on Mpc 
scales. Such correlations with radio sources are well docu- 
mented, z > 0.5 1 Jy Catalogue sources have been cross 
correlated with almost every available catalogue of 'low red- 
shift' extragalactic objects: Lick galaxies (Fugmann 1990, 
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Bartelmann & Schneider 1993) , IRAS galaxies, and EMSS 
(Bartelmann & Schneider 1994), ROSAT All-Sky Survey 
(Bartelmann, Schneider, & Hasinger 1994), Zwicky clusters 
(Seitz & Schneider 1995), APM galaxies (Bem'tez, N., & 
Martmez-Gonzales 1995, 1997), and Abell clusters (Wu & 
Han 1995). In all instances positive correlations were found, 
though at varying statistical significance, from 1.5 to >3 a. 
The large amplitude of the effect remained unexplained until 
recently. Dolag & Bartelmann (1997) and Sanz et al. (1997) 
reproduced the correlations on ~ 10' scales by incorporat- 
ing the non-linear growth of the matter power spectrum in 
the Universe. Including these mass fluctuations leads to an 
order of magnitude increase in the QSO-galaxy correlations 
on Mpc scales. 

On yet larger scales, 20' — 1°, there are three exist- 
ing studies with homogeneous complete catalogues of ob- 
jects. Rodrigues- Williams & Hogan (1994) looked at correla- 
tions between LBQS QSOs (Hewett et al. 1995) and Zwicky 
galaxy clusters, and found a significant signal with a subset 
of QSOs at 1.4 < z < 2.2. Seitz & Schneider (1995) extended 
the study to 1 Jy radio sources, and a range of QSO sub- 
samples based on redshift and apparent magnitude; sources 
at z ~ 1 were found to be associated with foreground clus- 
ters at 97.7% significance level. While Rodrigues- Williams 
& Hogan argue that the correlations could not be the result 
of patchy obscuration by the Galactic dust, the study of 
Seitz & Schneider, who used radio selected sources, proved 
so beyond a doubt. The physical nature of associations can 
be definitely ruled out based on the discrepant redshifts of 
QSOs and clusters, leaving weak gravitational lensing as the 
only plausible explanation. Ferreras et al. (1997) detect a 
strong QSO-galaxy anticorrelation between faint, z < 1.6, 
optically selected QSOs in a 5.5 deg 2 region close to the 
North Galactic Pole. The authors attribute the effect to the 
selected biases associated with identifying QSOs in crowded 
areas. 

In this paper we extend the previous work by examining 
the correlations between optically selected QSOs and fore- 
ground galaxies, on scales of up to a degree. In the following 
sections we measure the QSO-galaxy cross-correlation signal 
(Section^)), compare it to theoretical predictions (Section^), 
and discuss a few areas of cosmology that our observations 
will have a considerable impact on (Sections [|, and 0) . 



2 DATA 

The Large Bright Quasar Survey (LBQS, Hewett et al. 
1995) is the largest homogeneous catalogue of optically 
selected QSOs. Candidates were obtained from machine- 
scanned direct and objective-prism UK Schmidt Telescope 
plates, based on several selection criteria, including blue 
colour excess and presence of strong emission lines. The final 
list of QSOs was compiled after spectroscopic follow-up of 
the candidates at the Multiple Mirror and other telescopes. 
The LBQS contains over a 1000 QSOs between tub = 16 
and w 18.5 in 18 high Galactic latitude fields. The QSO 
redshift distribution, presented in Figure 9 of Hewett et al. 
(1995), is smooth and contains no gaps in the redshift range 
between 0.2 to 3.5. In the present work we only use the 
eleven equatorial LBQS fields. We do not use four fields in 
the direction of the Virgo Cluster, because the faint galaxy 
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Figure 1. Redshift distribution of the objects used in the present 
work. The three histograms are for QSOs with mg < 17.5, 18.0, 
and 18.5, all taken from eleven LBQS equatorial fields. The solid 
line is the estimated redshift distribution (Maddox et al. 1996) 
for 18.5 < mji < 20.0 APM galaxies (arbitrary vertical normal- 
ization). 



counts in the these fields can be severely contaminated by 
the Virgo galaxies. We do not use two fields in the direc- 
tion of the Galactic bulge, because of star contamination. 
We also did not use one field near the South Galactic pole. 

The APM Catalogue (Irwin et al. 1994) was compiled 
from the scans of Schmidt plates, carried out by the Au- 
tomatic Plate Measuring facility at Cambridge. The Cata- 
logue is a list of objects, detected on red and blue plates 
separately, and classified based on their morphology as star- 
like, extended, noise, or blend. The galaxy data we use is 
only approximately magnitude calibrated; two fields can be 
offset by as much as one magnitude. However, we do not at- 
tempt to correct for this because the difference is not large 
for our purposes, and our analysis compensates for it. 

Each Schmidt plate is ~ 6° across, but in order to re- 
duce vignetting problems we only use objects located within 
2.7° of the plate centre. 

Figure j] shows the redshift distribution of LBQS QSOs 
and the APM galaxies. The QSOs (three histograms for 
tub < 17.5, 18.0, and 18.5) are those found in 11 equa- 
torial fields used in this work. The solid curve represents 
18.5 < ma < 20.0 galaxies, and is based on the redshift dis- 
tribution estimated by Maddox et al. (1996). Because the 
latter analysis applies to APM galaxies detected on the blue 
plates while we mostly use red plates, we assumed a uni- 
form colour transformation of B — R — 1.5, which is the 
average colour of the 18.5 < thr < 20.0 galaxies detected 
on both red and blue APM plates. It is also consistent with 
the B — R colour derived by Metcalfe et al. (1995). The peak 
of the galaxy redshift distribution is at z « 0.2, and virtually 
no galaxies lie beyond z = 0.7. 
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3 QSO-GALAXY CORRELATIONS 

In this Section we ask ourselves if the sky-projected distribu- 
tions of the APM galaxies and background LBQS QSOs are 
correlated. The amplitude of the signal, if any, is expected 
to be small, therefore it is important to take into account 
any possible biases that may affect the signal. In particu- 
lar, Schmidt plates suffer from radial sensitivity gradients, 
mainly caused by vignetting, which results in radially depen- 
dent object number density. These gradients are small, and 
can be different for stars vs. galaxies, and for faint vs. bright 
objects. Our cross-correlation analysis, which we describe in 
the next Section, takes these effects into account. 



3.1 Method — cross-correlation estimator 

The standard measure of the clustering of objects is the two- 
point correlation function. The commonly considered case, 
the auto-correlation function of galaxies, has been explored 
in the literature in great depth; Peebles (1980), Hewett 
(1982), Hamilton (1993) have derived estimators that can 
tackle various observational limitations of the data, like fi- 
nite field size, average density variations on large angular 
scales, plate inhomogeneities, etc. Our case, QSO-galaxy 
cross-correlation, is different from galaxy-galaxy correlation 
not only because two different sets of objects are being con- 
sidered but also because QSOs, as opposed to galaxies, are 
rare; a typical field contains of the order of 10 4 galaxies down 
to niR — 20, but only a handful of bright QSOs. Therefore 
the various refinements on the standard estimator discussed 
in the literature are not applicable in the present case. 

The estimator that we use to compute the QSO-galaxy 
cross-correlation function is defined by, 



u(0) 



DqDg 
(R'qDg) 



1, 



(1) 



where DqDq is the actual number of QSO-galaxy pairs of 
a given separation, and R'qDg is the number of random 
QSO-real galaxy pairs. Random QSO positions are obtained 
for every real QSO by randomly 'scattering' it in the az- 
imuthal direction with respect to the field centre, while keep- 
ing its radial distance the same. This estimator deals suc- 
cessfully with the problem of low QSO numbers, our circular 
field boundaries, and radial sensitivity gradients on Schmidt 
plates. Had we used random QSOs that were scattered in 
both radial and azimuthal directions (or, equivalently in x 
and y), the presence of plate edges together with the small 
number of QSOs would have resulted in spurious correla- 
tions. Since azimuthally scattered random QSOs sample the 
same radius-dependent galaxy density on the plates as the 
corresponding real QSOs, radial plate gradients cancel out. 
To remind ourselves that our estimator is not the standard 
one, we use R'q, instead of Rq, in eq. |l|. The angular brackets 
in the denominator of eq. Q mean that we take the average 
of 100 random realizations for every real QSO. 

Note that for QSOs located close to plate centres, 
eq. [j] can underestimate the absolute amplitude of the cross- 
correlation signal. Suppose the plate centre happens to have 
a real excess or deficit of galaxies. All the randomly gen- 
erated QSOs close to field centre will also be sampling the 
same real excess or deficit of galaxies, and thus the cross- 
correlation signal will not be detected. However this is a 



small effect because only a small percentage of the QSOs 
are close to the plate centre. 

The cross-correlation as a function of angular scale, 6, 
is first calculated for every field separately. The final u(6) 
is the average of individual field contributions, weighted by 
the number of QSOs they contain. We do not attempt to 
estimate correlations on scales larger than a single plate, i.e. 
6>4°. 

To test our estimator we calculate the cross-correlation 
signal between galaxies (extended objects) and Galactic 
stars (star-like objects), both taken from the red APM 
plates. In this particular example galaxies are confined to 
a magnitude range between 18.5 and 20.0, while stars have 
magnitudes between 17.0 and 18.0. The surface density of 
these objects as a function of the distance from plate cen- 
tre, averaged over 11 plates, is shown in the top panel of 
Figure |i[ The empty circles are galaxies, filled circles are 
stars. Both types of objects show radial gradients. The star- 
galaxy cross-correlation estimated using eq. |l| is represented 
as star symbols in the lower panel of the same Figure. The 
cross-correlation estimated using the standard method, i.e. 
by scattering random QSOs in the x and y directions, is 
shown as square symbols. The estimator defined by eq. [t] 
has successfully compensated for galaxy and star radial gra- 
dients; its application shows that there is no correlation be- 
tween stars and galaxies, as expected. On the other hand, 
the standard estimator's results are dominated by spurious 
effects. We therefore use the estimator defined by eq. |l| in 
the rest of this paper. 

3.2 Results of cross-correlation 

Now we estimate QSO-galaxy correlations using one subset 
of QSOs, and different galaxy subsamples defined by their 
apparent magnitudes. As a check, for every QSO-galaxy cor- 
relation we calculate QSO-star correlation, with stars in the 
same magnitude range as galaxies. 

Figure |^ shows QSO-galaxy (top panel) and QSO-star 
(bottom panel) cross-correlation for QSOs with z > 1 and 
m < 18, and four galaxy subsamples, all taken from the red 
plates, with magnitude ranges: 
16.00 -* 19.50 (empty circles); 
17.32 -» 19.68 (triangles); 
18.50 -» 20.00 (solid circles); 
19.54 -> 20.46 (stars). 

The bin widths were chosen such that the number of galax- 
ies per bin is approximately the same in all bins. The QSO- 
galaxy cross-correlation signal, luqg{0), is detected clearly, 
while Galactic stars do not correlate with QSOs, as ex- 
pected. The signal persists to separations of ~ 75'; beyond 
that the slope of ujqg(6) becomes nearly zero. At #>150' the 
signal is mostly negative because of the integral constraint 
imposed on the correlation function. Though all four galaxy 
magnitude bins show comparable signal, the strongest cor- 
relations are with 18.5 < vtir < 20.0, i.e. the "faintest-but- 
one" bin. Perhaps this is not surprising: if the signal is due 
to weak lensing then fainter galaxy samples, i.e. those at 
higher redshifts, are expected to be better lenses for QSOs 
at z > 1. On the other hand, the faintest magnitude bin is 
probably contaminated by stars, which would tend to dilute 
the observed ujqg- 

Figure ^| is the same as Figure only the galaxies here 
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Figure 2. Top panel: normalized radial number density of 18.5 < 
mji < 20.0 APM galaxies around plate centres, summed over 11 
equatorial fields (empty circles), and 17.0 < m^j < 18.0 APM 
stars (filled circles). Both types of objects show radial density gra- 
dients. Bottom panel: cross-correlation between stars and galaxies 
using eq. [l] (star symbols) and 'standard' (triangles) estimators. 
The l-a errorbars, estimated as the standard deviation of the 
mean of 11 LBQS fields, are plotted in the latter case. Note the 
horizontal scale in the bottom panel is twice as long as in the 
top panel, since separations as large as the plate diameter can be 
considered. 



were taken from the blue plates. The magnitude ranges are: 
17.50 -► 21.00 (empty circles); 
18.82 -> 21.18 (triangles); 
20.00 -> 21.50 (solid circles); 
21.04 -> 21.96 (stars). 

The apparent association between 'stars' on blue APM 
plates and QSOs is due to a population of compact red 
galaxies, which appear extended on red plates, but are stel- 
lar on blue plates. This population has an average colour of 
B — R — 2.1, much redder than typical galaxies in the 18.5 
to 20.0 magnitude range detected on the red plates. The ex- 
istence of this population of galaxies is known, and prelim- 
inary spectroscopic work indicates that these are ellipticals 
at a typical redshift of 0.3. (Hewett et al. 1997). Recently, 
a faint lensed arclet of a high-redshift source has been ob- 
served centred on one of these galaxies (Hewett et al. 1997) . 
This population of compact red ellipticals at z ~ 0.3 is also 
an important contributor to the weak lensing induced cor- 
relation between QSOs and galaxies. Because these galaxies 
appear stellar on blue APM plates, a positive correlation be- 
tween QSOs and 'stars' is detected on these plates (Fig. 
bottom panel). Also, because this population is not detected 
as 'galaxies' on blue plates, the corresponding QSO-galaxy 
correlation (Fig. Q top panel) is weaker compared to that 
on the red plates. Because of these factors affecting the blue 
APM plates, we will only use red plates in the rest of the 
analysis. 

Before we proceed, it is important to note that the am- 



Figure 3. QSO-galaxy (top) and QSO-star (bottom) correla- 
tions using the same QSO subsample (z > 1, m < 18.0), and 
four galaxy subsamples taken from the red APM plates, with 
magnitude ranges: 16.00 — > 19.50 (empty circles); 17.32 — > 19.68 
(triangles); 18.50 -» 20.00 (solid circles): 19.54 -> 20.46 (stars). 
Error-bars are la standard deviation of the mean of 11 LBQS 
fields, and are shown only for the 18.50 < mji < 20.00 case. 
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Figure 4. QSO-galaxy (top) and QSO-star (bottom) correla- 
tions using the same QSO subsample (z > 1, m < 18.0), and 
four galaxy subsamples taken from the blue APM plates, with 
magnitude ranges: 17.50 — > 21.00 (empty circles); 18.82 — > 21.18 
(triangles); 20.00 -* 21.50 (solid circles); 21.04 -> 21.96 (stars). 
Error-bars are la standard deviation of the mean of 11 LBQS 
fields, and are shown only for the 20.00 < mji < 21.50 case. 
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Figure 5. Contour plots of cumulative ujqg and ujqs for 6 < 
75' (left), and the corresponding statistical significance (right). 
The correlation contours are drawn at 0.98, 0.99 (thick and think 
dashed lines), and 1.02, 1.01 (thick and thin solid lines) levels. 
The significance contours are at 90, 99% (thin and thick lines). 
Note that APM stars show no correlations with QSOs. 



plitude of the signal that we measure here is a lower limit to 
the true signal bec ause of three effects. First, as explained 
above (Section 3.1) our cross-correlation estimator can bias 
the amplitude of the signal low. Second, because QSO can- 
didates for the LBQS were selected using prism plates, can- 
didates in crowded areas, e.g. high galaxy density, are more 
likely to be rejected as their spectra have a higher chance of 
being 'corrupted' by a superimposed galaxy image. Third, 
APM object classifier and estimated object magnitudes are 
not perfect, and hence any true signal which is due to galax- 
ies in any given magnitude range will be diluted by stars and 
by objects of other apparent magnitudes. 



3.3 Integrated correlations for 8 < 75' 

In the previous Section we used one subsample of QSOs. To 
eliminate any bias related to that particular choice of QSO 
properties, we extend the QSO-galaxy cross-correlation 
analysis to a range of QSO subsamples each having a mini- 
mum redshift and limiting apparent magnitude, ZQ, m in and 
TUQ^im- For each subsample we calculate the integrated cor- 
relation amplitude within 6 — 75', i.e. we characterize each 
{zQ, m in, rnQ,u m ) subsample by one number. This integrated 
correlation amplitude is plotted as a contour plot in Figure ^| 
The top two panels of Figure |E| show the contours of 
constant QSO-galaxy correlation amplitude (left panel) and 
statistical significance (right panel); while the bottom two 
panels are the same, but for QSO-star correlations. The cor- 
relation contours are drawn at 0.98, 0.99 (thick and thin 
dashed lines), and 1.02, 1.01 (thick and thin solid lines) lev- 
els. The significance contours are at 90 and 99% (thin and 
thick lines). The significance level is the fraction of synthetic 



QSO subsamples, out a total of 100, that show correlations 
weaker that the corresponding real QSO subsample. This 
estimate makes no a priori assumptions about the distribu- 
tion of errors. The right corner of each panel is shaded where 
the number of QSOs per subsample is less than 10, so even if 
the correlations appear significant, they should be regarded 
with caution because not enough QSOs are being averaged 
over. 

It is apparent from Figure ^ that QSO-galaxy correla- 
tions are statistically significant for some QSO subsamples, 
whereas no significant correlations, above 90% c.L, are de- 
tected between QSOs and Galactic stars. Because of the null 
result for QSO-star and star-galaxy (Figure ^) correlations, 
we conclude that QSO-galaxy signal is real, and not an ar- 
tifact of plate sensitivity gradients, etc. 

What are the possible causes of these correlations? The 
correlations cannot be due to physical QSO-galaxy associ- 
ations, because of the QSO/galaxy redshift mismatch. An- 
other possibility is that patchy Galactic obscuration 'creates' 
galaxy and QSO over-densities, which are then necessarily 
correlated. This hypothesis can be ruled out because the 
observed correlations show a strong variation with ZQ, m i n 
both in amplitude and significance, whereas dust would not 
be able to differentiate between QSOs at different redshifts. 
Furthermore, the fluctuations in the distribution of Galactic 
dust are not extreme enough to produce the observed varia- 
tions in galaxy and QSO number densities. The rms fluctua- 
tions in obscuration in these high Galactic latitude (|6|>45°) 
LBQS fields is typically 0.1 magnitudes in B (Burstein & 
Heiles 1978, Schlegel et al. 1997) which would produce rms 
fluctuations in projected galaxy density of a factor of 1.07, 
while the observed rms value is 4 on ~ 10' scales, and 1.4 
on degree scales. Thus, the observed galaxy density fluctua- 
tions in these APM fields are primarily due to galaxy clus- 
tering, and not patchy Galactic dust obscuration. Therefore 
we conclude that the correlations are not due to dust. 

The only remaining explanation is the weak gravita- 
tional lensing of the background QSOs by the matter asso- 
ciated with the APM galaxies. This interpretation is sup- 
ported by the trends seen in the top panels of Figure [| The 
correlations are strongest with QSOs at z<;l, consistent with 
theoretical expectations (Dolag & Bartelmann 1997, Sanz et 
al. 1997). The apparent magnitude-dependent behavior of 
the correlations is also consistent with the lensing interpre- 
tation; ujqg is strongest for m,Q<18. At fainter magnitudes 
the slope of the QSO luminosity function becomes shallower 
and so amplification bias, and hence correlation strength, 
become less pronounced. At magnitudes much brighter than 
~ 18 the significance of the QSO-galaxy correlations is low 
due to small QSO numbers. 

As an example of the actual ujqg as a function of sepa- 
ration, Figure ^| shows the correlation between z > 1 QSOs 
and galaxies, for three different mg.iim: 17.5, 18.0 and 18.5. 
The star symbols are QSO-star, and solid dots are QSO- 
galaxy correlations. The error-bars are la standard devia- 
tion of the mean of eleven LBQS fields. The results shown 
in the middle panel, itiq < 18, are the same as solid circles 
in Figure ^[ 



© 1998 RAS, MNRAS 000, 000-000 



6 L. L. R. Williams and M. Irwin 




Figure 6. QSO-galaxy (solid dots and lines) and QSO-star (star 
symbols) cross-correlation function. The galaxies are taken from 
the red APM plates, and have magnitudes between 18.5 and 20.0. 
Results for three QSO subsamples are presented, as indicated in 
the upper right corner of each panel. Note that the amplitude of 
u)qg increases rapidly for brighter QSOs. 



4 ANALYSIS AND COMPARISON WITH 
PREDICTIONS 

Having measured the weak lensing induced two-point corre- 
lation function between QSOs and galaxies, we now compare 
it to the oreti cal expectations, using two approaches (Sec- 
and 4.2) that differ mostly in the method used to 



tions 4.1 



estimate the dumpiness of the mass responsible for lensing. 
The results of these two approaches are compared to each 
other and to the observations. As observations, we chose the 
correlations presented in the middle panel of Figure [], i.e. 
for a QSO subsample with zq > 1 and iuq < 18.0. We plot 
these as filled circles in Figure]^, for separations where the 
amplitude of loqg remains consistently positive. 



4.1 Phenomenological Predictions 

The main assumption here is that the mass distribution re- 
sponsible for lensing is traced by the observed galaxies in the 
nearby Universe. The total 'thickness' of the lensing material 
can be estimated from the galaxy redshift distribution and 
an assumed Qrj> while the distribution of fluctuations on any 
given scale are obtained from the counts-in-cells distribution 
of the APM galaxies. 

First, we estimate the amount of matter associated with 
the observed APM galaxies, i.e. we calculate the optical 
depth for sources located at the redshift of the QSOs, 



TAPM = Pcrit^O 



(cdt/dz)(l + z) 

^crit (-2, Z s ) 



■dz. 



(2) 



Here, cdt is the thickness of the lensing slice at redshift z, 
p C rit&o{l + z) 3 is its mass density, and T, cr it{z, z s ) is the crit- 
ical lensing surface mass density at z for a source at z„. The 



Figure 7. Observed and predicted QSO-galaxy correlations; note 
the log-log scale. The different lines and points are described in 
the Figure. The solid do ts re pres ent the z > 1 and m < 18.0 QSO 
subsample. See Sections 4.1 and 



for details. 



latter is assumed to be 1.5, the average redshift of our QSO 
subsample. We assume an open Universe with flo — 0.3. The 
upper limit of integration, z max is taken to be the redshift at 
which the APM galaxy redshift distribution (Figure |l]) drops 
to half of its peak value; z max is about 0.33. In other words, 
we assume that the optical depth of the mass traced by the 
APM galaxies is the total optical depth up to a redshift of 
0.33. Since the particular choice of z max is somewhat arbi- 
trary, we also quote results for z max =Q.2& and 0.39, corre- 
sponding to redshifts where the APM galaxy redshift distri- 
bution drops to 3/4 and 1/4 of its peak value, respectively. 
The corresponding range of tapm values is from 0.012 to 
0.021, while for z max =0.33, it is 0.016. 

Next we assume that on scales larger than a few ar- 
cminutes the APM galaxies trace the total matter distri- 
bution (up to z max ) with a biasing factor b. The projected 
galaxy clustering is described by the scale-dependent counts- 
in-cells distribution of the APM galaxies, p(u\ 6) da, where 
a is the galaxy number density normalized by the average 
density. A patch of sky with density a produces an ampli- 
fication A{o) « 1 + 2tapm{p — l)/&, in the weak lensing 
limit. The corresponding over-density of QSOs background 
to this patch is given by, 



q(a) = A 



1 + (2tapm[(t -l]/b) (2.5a - 1) 



(3) 



where a is the slope of the QSO number counts; (a) = 1.1 
for our QSO subsample. 

The cross-correlation between QSOs and galaxies is 
then estimated as 



(4) 



wqg(< + l -/ p{a\0)q{a)ad(J. 



The galaxy auto-correlation function is given by an 
equation similar to eq. M, 
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logg(<6) + 1~ p(a\6)a 2 da. (5) 
Jo 

The ratio of ujqg to ujgg can be derived analytically, in the 
weak lensing regime. Making use of the fact that p(cr\ 9) is 
a normalized probability distribution with an average value 
of a equal to 1, the ratio of eq. ^| and eq. ^ gives, 

ujqg(6) « (2r A p M /b) (2.5a ~ 1)uj gg (0), (6) 

This equation applies if the same galaxy counts-in-cells dis- 
tribution is used to calculate ujqg(0) and UJ GG (0), and gives 
the QSO-galaxy correlation amplitude due to lensing by 
these galaxies only. 

For tapm = 0.016 (i.e. z max = 0.33), b = 1, 
and (a)=l.l, the ratio ujqg/uJgg ~ 0.056. Galaxy auto- 
correlation function from the APM galaxies was calculated 
by Maddox et al. (1990), and is plotted as dot-dashed line 
in Figure [?]. The amplitude of ujgg is appropriate for the 
galaxy magnitude range used here. With this ujgg, QSO- 
galaxy cross-correlation function should lie 1.25 (in the log) 
below ujgg- Since the observed ujqg is roughly of the same 
amplitude as the galaxy-galaxy correlations, the predictions 
based on eq. ^| are a factor of 15 below observations. 

As a check, eqs. § and | can be also applied directly 
to the counts-in-cells data to calculate the correlation func- 
tions at a range of separations, assuming the optical depth 
of the APM galaxies. The results are shown in Figure m, as 
empty squares (ujgg) and empty circles (ujqg)- The vertical 
offset of the two is about a factor of 15, in agreement with 
eq. |(| Because vertical normalization of ujgg scales with the 
magnitude of the galaxies, and the data we use was not 
properly magnitude calibrated, we have adjusted the verti- 
cal normalization of these two estimated functions such that 
our derived ujgg matches that of Maddox et al. (1990). 

The empty circles assume z max =0.33; the dotted line 
through them is a rough fit. The upper and lower dotted 
lines are for z maa; =0.39 and 0.28 respectively, and bracket 
the range of our predicted values. The model predictions 
cover an angular scale range of 4' to 1°; on scales smaller 
than 4' the number of galaxies per cell is small, and so the 
resultant counts-in-cells distribution is dominated by Pois- 
son noise. Scales much larger than 1° are comparable in size 
to the Schmidt plates, so the derived counts-in-cells will be 
substantially narrower then the true distribution, because 
we normalize the galaxy number density to the plate's av- 
erage, on each plate separately, thus ignoring any power on 
larger scales. For the same reason, the slope of our ujgg(&) 
is somewhat steeper than the commonly derived value of 
7 w 0.7 (Maddox et al. 1990). 

As mentioned above, our phenomenological model pre- 
dictions fall short of the observations by about a factor of 
10-20. Is our model too simplistic? 

4.2 Theoretical Predictions from the Literature 

Dolag & Bartelmann (1997) and Sanz et al. (1997) used 
an alternative way to derive mass fluctuations in the Uni- 
verse. Starting with a particular type of cosmology and an 
initial matter power spectrum, they used the results of a 
semi-analytic prescription for the non-linear evolution of the 
power spectrum (Peacock & Dodds 1996) to derive the spec- 
trum of mass distribution as a function of redshift and scale. 



Assuming weak lensing regime, slope of QSO counts a, and 
a biasing parameter, the mass distribution spectrum was 
used to predict the observed QSO-galaxy correlations. These 
calculations were primarily motivated by an observed corre- 
lation between PKS radio-selected QSOs and APM galax- 
ies, as derived by Bem'tez & Martmez-Gonzalez (1997), and 
shown as star symbols in Figure ^. Within their range of 
applicability, i.e. 6 > 1 — 2', both analytical models re- 
produce the PKS- APM results well. Correlations on scales 
larger than 15' could not have been detected by Benftez & 
Martmez-Gonzalez because they restricted their analysis to 
patches of 15' radius around each QSO. 

To compare these two model predictions to our observed 
ujqg(&), we rescaled their results to a = 1.1. Dolag & Bartel- 
mann results were additionally rescaled from h — 0.7 to 
h = 0.5 using their Figure 4. Thepredictions are plotted 
as dashed and solid lines in Figure fit both are for Q, = 0.3, 
A = 0. Similar to the results of the phenomenological model, 
the analytical models underestimate the amplitude of ob- 
served ujqg by a factor of 10 on large angular scales. 

Since phenomenological and analytical models used dif- 
ferent routes to arrive at ujqg, it is encouraging that they 
agree reasonably well with each other. A factor of 1.5-2 dis- 
crepancy (solid and dashed lines versus dotted lines) arises 
probably because our phenomenological model ignored lens- 
ing by structures at z>0.35, which are not sampled by galax- 
ies in the APM Catalogue, but are still efficient lenses for the 
QSOs. Turning the argument around, the small discrepancy 
means that nearby galaxies, z<0.35, dominate the weak lens- 
ing of z > 1 QSOs. 



5 DISCUSSION 

In the last Section we have seen that the amplitude of ob- 
served QSO-galaxy correlation is under-predicted by theo- 
retical models. The discrepancy is a factor of 3-4 at 8^10' , 
and increases to a factor of 10 on degree scales. 

Can a reasonable change of parameters account for such 
a large discrepancy? Assuming that the effects of weak lens- 
ing are dealt with correctly, and the signal is real and is due 
to lensing, we are left with two possibilities (see eq. ^|); either 
the slope of the QSO number counts at z > 1 is very steep, 
or mass fluctuations on relevant scales, i.e. between a few 
and 15 Mpc, are more extreme than fluctuations in galaxy 
number density. The amplitude of ujqg is quite sensitive to 
the slope of the QSO number counts; however to reproduce 
the observations, a would have to be ~ 8, while the ob- 
served slope values for LBQS QSO in this redshift range 
lie between 0.8 and 1.5, comparable to the QSO luminosity 
function slopes, ~ 1.6, derived from the LBQS (Hewett et 
al. 1993) in this redshift range. 

Alternatively, we could require that the galaxies are 
anti-biased with respect to the mass by a factor of ~ 10, 
(b ~ 0.1) implying a as value of around 10, which is hard 
to reconcile with any gravitational clustering scenario, and 
observations of the large scale structure. The value of fl 06 as 
is well constrained by the abundance of galaxy clusters; 
White, Efstathiou & Frenk (1993) estimate Q 0f V 8 to be 
~ 0.6. Based on cluster peculiar velocities, Watkins (1997) 
derives f2 a6 <T S to be 0.44. Thus an as ~ 10 would re- 
quire an Q so low, it would just agree with tlbaryon from 
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primordial nucleosynthesis constraints (Hogan 1996), but 
would strongly contradict dynamical measurements on clus- 
ter scales (£7=0.24, Carlberg et al. 1996), and power spec- 
trum shape parameter constraints (r=Q /i=0.2-0.3, Mad- 
dox et al. 1996). Furthermore, such a high a$ value is ruled 
out by the recent direct estimation by Fan, Bahcall & Cen 
(1997), who compute as = 0.83 ± 0.15 from the redshift 
evolution of galaxy clusters. 

A combination of factors, for example, 6=0.3 and a=2.3 
shares the burden between the two observables, but is still 
quite unpalatable. In fact, it is difficult to see how any rea- 
sonable change of parameters or assumptions can bring these 
up to the observed amplitude of loqg on ~ degree scales. 

If the explanation does not lie with either the lensers 
or the lensees, then maybe the lensing process has to be 
looked at more carefully. It is conceivable, for example, that 
some hitherto unexplored non-linearities in the description 
of light propagation through a clumpy universe could be 
responsible. 

At present, the factor of 10 discrepancy between obser- 
vations and predictions remains unaccounted for. 



6 IMPLICATIONS FOR OBSERVATIONAL 
COSMOLOGY 

The results of the correlations described above have interest- 
ing implications for the observed distribution and properties 
of certain types of high-redshift sources. Here we explore a 
few of these. 

The results imply that the projected density of fore- 
ground structure affects the density of high redshift objects. 
To quantify this statement in the case of QSOs our results 
can be represented as a plot of QSO number density as a 
function of foreground galaxy density on a given angular 
scale. To do that, we randomly lay down a large number of 
circular patches over LBQS fields. In each patch we calcu- 
late galaxy number density. For all patches of a given galaxy 
density we then calculate the average QSO density, and re- 
peat the process for several values of galaxy density. The 
results are shown in Figure ^, using patches of radius 12', 
and three QSO subsamples, all with zq > 1. The error-bars 
are la standard deviations of the mean of 10 different real- 
izations of the experiment just described. 

It is seen that the observed projected density of QSOs 
is a function of the foreground galaxy density, in the sense 
that bright QSOs are found preferentially in the directions of 
over-dense galaxy regions; for example, a tuq < 17.5, zq > 1 
optically selected QSO is twice as likely to be found in the 
direction where the galaxy density is 1.5 times the average. 
However, QSOs are not the only class of objects affected 
by weak lensing on large scales. Any population of high- 
redshift objects, whose number counts slope is substantially 
different from 0.4 will be affected in a similar fashion. Intrin- 
sically bright galaxies, those on the exponential part of the 
Schechter luminosity function (Schechter 1976), confined to 
a narrow redshift range, are an example. As indicated in Fig- 
ure ^, both under- and over-densities of background sources 
can be expected. Thus, this bears direct relevance to the 
Hubble Deep Field (R. E. Williams et al. 1994), which was 
intentionally chosen to lie in the direction of the sky devoid 
of nearby, z<0.3, structure (the nearest galaxy cluster is 48' 



away, etc). Depending on the exact under-density of nearby 
galaxies in the HDF, and the intrinsic luminosity of high- 
redshift galaxies, the latter can be depleted by a factor of 
up to 5. 

A further implication is that the rms dispersion in 
the observed luminosities of standard candles, and observed 
sizes of standard rulers at high-redshift would be somewhat 
larger than currently predicted from numerical simulations 
of the lensing effects (Wambsganss et al. 1997). A rough ex- 
trapolation from our results shows that the rms dispersion 
in amplifications would be 0.1 — 0.2 magnitudes, comparable 
to the light-curve shape corrected dispersion in the luminosi- 
ties of four z ~ 1 supernovae of Type la (Garnavich et al. 
1997), but smaller then the dispersion in sizes of compact ra- 
dio sources (Kellerman 1993) and double-lobed radio sources 
(Buchalter et al. 1997). The fact that the rms dispersion in 
these observed quantities is comparable to or larger than 
what is implied by our weak lensing analysis, means that 
the rms spread in lensing induced amplifications is not use- 
fully constrained by the current samples of standard candles 
and rulers. Because the weak lensing induced dispersion in 
the observed quantities is not Gaussian (or even symmetric), 
the overall trend with redshift can also be affected (also see 
Wambsganss et al. 1997). 

Another aspect of extragalactic astronomy that is af- 
fected by weak lensing is the density and redshift evolution 
of Ly-a forest clouds. It is evident from Figure ^Jthat bright, 
high redshift QSOs do not sample random lines of sight, as is 
generally hoped in the studies of intervening QSO Ly-a ab- 
sorption lines. Since the most efficient lenses for a wide range 
of source redshifts are at Zj<0.6, one would expect a corre- 
sponding increase in the detected density of high equivalent 
width Ly-a forest clouds, which are known to be loosely 
associated with galaxies (Stocke et al. 1995, Bowen et al. 
1996, Lanzetta et al. 1995) at these redshifts. An increase in 
dN /dz at z<l compared to a power law extrapolation from 
higher redshifts has been detected (Bahcall et al. 1993, Im- 
pey et al. 1996); a substantial contribution to this increase 
is probably due to weak lensing of back-lighting QSOs. 



7 CONCLUSIONS 

We have found that bright, high-redshift, optically selected 
LBQS QSOs are positively correlated with foreground APM 
galaxies. The cross-correlations are significant and are not 
physical in nature, or due to patchy Galactic obscuration. 
The most plausible explanation, which is also supported by 
the qualitative behavior of the correlations, is weak gravi- 
tational lensing of the background QSOs by the intervening 
dark matter traced by the APM galaxies. However, the am- 
plitude of correlations on degree scales, or ~ 10-15 Mpc at 
the redshift of the lenses, is a factor of 10 higher than pre- 
dicted from models. The discrepancy is hard to reconcile 
with what we currently believe to be true about the Uni- 
verse. 

The implications of the effects of weak lensing for the 
observed high redshift Universe are far reaching, therefore it 
is imperative to study such correlations further using other 
large uniform data sets. An analysis with radio selected high- 
redshift sources would be ideal, as these sources are immune 
to the effects of Galactic dust obscuration. Samples of high- 



© 1998 RAS, MNRAS 000, 000-000 



Weak tensing: LBQS and APM 9 




normalized galaxy density 



Figure 8. Normalized QSO density for a given normalized galaxy 
density, in circular patches of radius 12'. See Section ^ for details. 
Galaxies are 18.5 < mjj < 20.0; QSOs are at Z. > 1. Bright 
QSOs are more likely to be found in the directions over-dense in 
nearby galaxies, and avoid direction of lower than average galaxy 
densities. 

redshift sources, both optical and radio, observed down to 
faint flux limits would be useful, as both positive correlations 
and anti-correlations with foreground galaxies are expected, 
depending on the limiting flux of the source subsample. 

Finally, it is interesting to note that the QSO-galaxy 
correlations described here are probing the physical scale 
where the galaxy power spectrum is observed to have a 
kink, corresponding to a primordial feature in the true linear 
power spectrum, which is not reproduced by any variants of 
the CDM model (Peacock 1997). 
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